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Abstract 


It is well known that the resonator geometry strongly influences the resonant frequencies 
of an acoustical resonator and the generated nonlinear standing pressure waveform. Maximizing 
the ratio of maximum to minimum gas pressure at an end of an oscillating resonator by 
optimizing the cavity contour is investigated numerically. A quasi-Newton type scheme is used 
to find optimized axisymmetric resonator shapes to achieve the maximum pressure compression 
ratio. The acoustical field is solved using a one-dimensional model, and the resonance frequency 
shift and hysteresis effects are obtained through an automation scheme based on continuation 
methods. Results are presented from optimizing cone, horn-cone, and cosine resonator 
geometries. Significant performance improvement is found in the optimized shapes over others 
previously published. Different optimized shapes are found when starting with different initial 
guesses, indicating multiple local extrema. The numerical model is validated by comparing with 
the experimental results of a horn-cone shaped resonator. 

PACS numbers: 43.25.Gf, 43.25.Cb 


2 



Li, JASA 


I. INTRODUCTION 

The waveform of the standing wave in an oscillating closed cavity is strongly influenced 
by the geometry of the resonator cavity. It is well known that shocks form in a cylindrical tube 
when the interior gas is oscillating at its resonant frequency. Lawrenson et al . 1 at MacroSonix 
Corp first exploited the shape dependence and obtained high-amplitude and shock-free acoustic 
pressures in axisymmetric tubes of varying cross sections, referred to as resonant macrosonic 
synthesis (RMS). Peak acoustic pressures that measure three to four times ambient pressure and 
maximum to minimum pressure ratios of 27 were observed in shaped cavities. The size of the 
demonstrated overpressure reached the level that is required by commercial applications such as 
acoustic pumps or compressors. The researchers considered these types of axisymmetric shapes: 
cylinder, cone, hone-cone, and bulb. They concluded that the hone-cone resonator shape 
generated the highest overpressure for a given input power. They also demonstrated that the 
overall characteristics of the waveform does not change when a resonator is filled with different 
gases. 

A companion paper by Ilinskii et al? developed a one-dimensional frequency-domain 
model for studying the RMS numerically. The results confirmed the nonlinear standing 
waveform and the related characteristics such as shape-induced resonance hardening and 
softening observed in the experiments by Lawrenson et al? To account for the energy losses in 
the boundary layer along cavity wall, Ilinskii et al 3 later modified the one-dimensional model by 
introducing an additional term in the continuity equation and used a turbulence model. Hamilton 
et al. 4 analytically investigated the relationship between the natural frequency of a nonlinear 
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acoustic resonator and its shape as well as the nonlinear interactions of modes in the resonator. 
Chun and Kim 5 numerically investigated cosine shaped resonators in addition to cylindrical and 
conical shapes using high-order finite-difference approximations. They concluded that the half 
cosine-shape is more suitable to induce high compression ratio than other shapes under certain 
assumptions. Recently, Erickson and Zinn 6 used the Galerkin method to solve the one- 
dimensional model and found a non-monotonic increase in compression ratio when the flare 
constant is raised for a class of hom-shaped resonators. To serve commercial needs, such as in 
acoustic gas compressor and acoustic liquid pump, the objective is to find an optimized shape for 
generating higher overpressure. However, the optimization procedure and results have not yet 
been discussed. 

In this article, the numerical schemes are introduced and the results are presented for 
optimizing the shape parameters that yield the highest maximum-to-minimum pressure ratio in 
each of the following resonator shapes: cone, hom-cone, and cosine-shape. In Sec. II, the 
modeling equations are presented; in Sec. Ill, the numerical schemes are described; in Sec. IV, 
the numerical model and experiment results are compared and the resulting optimized shapes are 
discussed. 


II. GOVERNING EQUATIONS 

In this section, the one-dimensional model for computing the acoustic wave field in an 
axisymmetric resonator is briefly described. The equations are presented for completeness and 
the details of the derivation are given in the work by Ilinksii et al. 2 . 
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Consider the acoustic field in an oscillating resonator of length / driven by an external 
force. The resonator is axisymmetric with the inner radius given by r=r(x), 0<x<l, where x is 
the coordinate along the axis of symmetry. 

The density of the gas p, the velocity u, and the pressure p satisfy the conservation of mass 


3e + 4i. (r » = o, 

dt r 1 dx K H ’ 


( 1 ) 


and the conservation of momentum 


ou ou 


1 dp 


— + u — = —~a(t) + 

dt dx p dx 


(£ + 477 / 3 ) d 
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r 2 dx K 


( 2 ) 


p dx 

where a(t) is the acceleration of the resonator enforced by the external force; ^ and r| are 
coefficients of viscosity, rj is the shear viscosity, and £, is the bulk viscosity that results from 
nonequilibrium deviations between the actual local pressure and the thermodynamic pressure 
The state equation is specified by that of an ideal gas 


P = Po{p / PoY > 


( 3 ) 


where p 0 is the ambient pressure, p 0 is the density, and y is the ratio of specific heats of the gas. 

The no-penetration boundary conditions at the two ends require that the velocity vanish at 
x = 0 and / . 


The quasi-one-dimensional compressible Navier-Stokes equations (l)-(3) can be solved 
using different numerical methods. Instead of finite difference method or finite element method, 
the frequency-domain method is employed, where the unknown variables are expressed in terms 
of finite Fourier series in time and the Fourier coefficients are solved for. Both the experimental 
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and numerical results show that the time harmonics of the dependent variables, such as the 
pressure p , decay rapidly as the frequency increases. This ensures that the number of time 
harmonics, N , needed for accurate results is small. In our numerical simulations, we have set 
the value of N to be ten and found that increasing N does not change our results. 

Following Ilinksii et al 2 by expressing the variables in Fourier series, the Eqs. (1) and (2) 
can be reduced to a system of Ordinary Differential Equations (ODEs) for the Fourier 
coefficients of the velocity potential tp, defined as u-W(p: 


d< PL-2t 


dx r 


2 ’ 


/ dx 


(4) 


a N 

where v = r 2 — , and <p k , v k and a k are the Fourier coefficients defined as <p= V <p k e' kl , 

dx 


k=-N 


/V 

v= X v k e ' k ‘ , a = ^ **k elk * • The detailed expressions for D kl = D kl (v k ,<p k ,x) and 


k=-N 
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where 5 - ^ + ^ T 1 (2 . j s the viscosity, c 0 = . r~^- is the reference speed of sound, D' m is given 
Po V Po 


by D' m = -{y-\)xa m -im{y-\)v m -^-[v 2 ] m , and [v 2 ) k = £v w v,. The no-slip boundary 


l=-N+k 
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conditions at the two ends are translated to the equivalent conditions in Fourier space: v k = 0 
at* = 0,1 . 


After the Fourier coefficients { cp k } and { v k } are obtained by solving the boundary-value 
problem Eq.(4), the velocity potential (p and the modified velocity v are computed from the 
inverse FFT. The density p is given by the momentum equation 

i 


P_ 

Po 


1 - 





5 dv 
+ ax — ^ — 
r 2 dx 


V- > 


-i/ 


and the pressure p can be obtained from the state equation (3). 


(5) 


III. NUMERICAL METHODS 

In this section, the procedures for finding the optimal shape parameters are described so 
that the desired pressure waveform can be obtained. The overall steps are as follows: first, given 
a resonator shape and a fixed value of input power, the resonant frequency of the resonator is 
searched and the compression ratio, defined as the maximum-to-minimum pressure ratio at the 
narrow end, is computed; second, the optimization step is performed, which yields the next 
candidate for the optimal resonator shape; then, the first and the second steps are repeated until 
an optimal design is reached. 

To simplify the discussion, the following dimensionless variables are introduced, 
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X = - , T = cot, 
l 





( 6 ) 


where co is the frequency of the periodic force acted on the resonator, o)o = nc 0 / / is the 
fundamental frequency of a cylindrical resonator of length / . 


In this work, the acceleration of the resonator induced by external force is assumed to be 
harmonic, A(T) = Acos(T) . The physical parameters that determine the acoustic wave field in the 


resonator are: the acceleration amplitude A , the ratio of specific heats y, the attenuation 
coefficient g = — » , and the resonator oscillating frequency q - XL . 


c oPo 


co n 


A. Shape optimization 

Suppose the resonator shape R(X) is determined by a number of shape parameters, 
S 0 ,S v ...,S n For example, a cone is given by R{X)~ S 0 + . The goal is to find the best shape 

producing the highest value of the pressure compression ratio R c , which is defined as the ratio of 
maximum pressure to the minimum pressure at the narrow end of the resonator, i.e., 

RciSoA sj = Z^, a tx=o. ( 7 ) 

P min 

The compression ratio R c is a function of the shape parameters, the dimensionless frequency Q , 
and the history of Cl (due to the existence of hysteresis effects). The method for obtaining R c for 
a fixed resonator shape is explained later in the section. 
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Now, the method for finding the optimal resonator shape is presented. First, the practical 
conditions and constraints in the optimization are specified. The original experiments by 
Lawrenson et al .' and some commercial applications use refrigerant R-134a as the gas inside 
resonators. We set the specific heat ratio y = 1.2 in order to match the gas properties of the 
refrigerant. The viscosity -related parameter G is fixed at 0.01, and will be shown later to be the 
value that makes the results from our numerical model match with those from the physical 
experiments. This same value of G was found appropriate in the study of Illinski 2 . For the 
current optimization, the dimensionless radius of the resonator at the narrow end, R(X = 0) , is 

restricted to be greater than 0.020833 so that real resonators could be built using these optimized 
dimensions. 

The objective of the optimization is to maximize the compression ratio that would be 
observed in a real oscillating acoustic resonator. The electrodynamic shaker system used in this 
study and others 1 ' 6 has limited acceleration capacity that is reduced by adding weight to the 
resonator. If we used the constraint to be constant acceleration amplitude, the optimization 
scheme would predict an optimum conical resonator shape that has a very small narrow end and 
very large wide end 7 , resulting in heavy weight experimental hardware. However, in an 
experimental setup, this resonator would not produce high compression ratios due to the low 
acceleration generated by the shaker. A greater compression ratio would be obtained from 
lighter weight hardware oscillated at higher acceleration amplitude. Therefore, the amount of 
input power that is delivered to the resonator oscillation is fixed during the optimization, 
effectively favoring lighter weight resonators and more closely modeling the real system. 
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For a fixed amount of input power W , , the acceleration amplitude A from the mass of 
the resonator is deduced through the formula in dimensional variables 

(W,) = (p 0 col)-MA\ (8) 

n 

where M is the dimensionless mass of the resonator. Assuming the width and the density of the 
resonator wall are constant, the mass of the resonator M is proportional to the total surface area 
of the resonator. Equation (8) can be derived by assuming that most of the power consumed, 
<W|>, is due to the inertial force F required to oscillate the resonator at the frequency o . In 
dimensional variables, 


/Uf\ 4 f’"* A, 4 fV 4 A, 4 fV 4 du A . 

(IV.) = — I Fu at = — I mau at = — m — u dt = 

1 T ’ * T •» T » dt 


\TJ4 


= — mu 


( 9 ) 


l (=0 


where T p is one period of the oscillation. Equation (8) follows as the dimensionless velocity of 
the resonator is given by A cos T and the power is nondimensionalized by p 0 eol . 

Second, a quasi-Newton method, BFGS (Broyden-Fletcher-Goldfarb-Shanno 8 ), is used 
for maximizing the multi-variable nonlinear function R c (S 0 ,S l S n ) . Since the evaluation of the 

objective function R c itself involves solving a nonlinear system of ODEs, Eq. (4), many times, 
the gradient information of R c required for the BFGS method is not available analytically and is 

derived by partial derivatives using a numerical differentiation method via finite differences. 

This entails perturbing each design variables, S„ in turn and calculating the rate of change in the 
objective function R C (S 0 , 5, S n ) . For two shape parameters, the optimization takes 4 to 48 hours 
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of CPU time (depending on the type of the resonator shape and initial guess for the shape 
parameters) on a 1 ,3GHz Athlon T-Bird PC with the Lahey-Fujitsu FORTRAN compiler. 

B. Boundary value problem 

For a given shape of the resonator, the boundary value problem Eq. (4) is solved 
numerically by a Multiple Shooting Method. Since the amplitude of the pressure in a resonator 
strongly depends on the oscillating frequency Q , the quantity to be optimized is the maximum 
pressure compression ratio R c over the entire range of Q for a given resonator shape. Because 

of the hysteresis effects, the solution is not unique near the resonant frequency and the Multiple 
Shooting Method will not converge unless a good initial guess of the solution is provided. To 
circumvent the difficulty, a continuation method is implemented: the system of ODEs is solved 
starting from a frequency Q that is far away from the resonant frequency and the solution is used 
as an initial guess for solving the ODEs for increased or decreased Q; the steps are repeated until 
all branches of the solution for all values of Q near the resonance is completed. The maximum 
ratio R c among different values of Q is chosen as the compression ratio for the resonator. 

IV. RESULTS 

A. Characteristics of the standing waves 

Before presenting the results of the optimal resonator's shapes, some of the important 
properties of the pressure wave in a non-cylindrical resonator are illustrated. 

The horn-cone resonator shape is described by: 
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R(X) = 


fS 0 cosh(S,A), 
\a + j3X, 


forO<X<S 2 

forS 2 <X<l 


( 10 ) 


where 5 0 = 0.028333,5, =5.7264, S 2 = 0 . 25 , a = S 0 cosh(5, 5,), and = S 0 S ] sinh(S,S 2 ). The shape 
parameters for the homcone are obtained from those used in the experiments by Lawrenson et 

al ] 


In Fig. 1 , the pressure waveforms at the ends of a cylindrical resonator and the homcone 
resonator are shown at their corresponding resonance frequencies for the effective viscosity 

G-0.01 , the gas specific heat ratio y— 1 .2, and the acceleration A = 1 x 1(T 3 . Throughout this 
paper, the effective viscosity G and the gas specific heat ratio are held constant. At the same 
acceleration magnitude, the waveform at the narrow end of the homcone shows large variation in 
which the ratio of the maximum and minimum pressure, R c , exceeds the value of 30.5. On the 
other hand, the waveform for the cylinder shows the formation of shocks at resonance and the 
compression ratio R c is below 1.32. The difference in the waveforms shows the strong 

dependence of the acoustic field on the geometry of a resonator. The graph also shows that the 
variation in pressure at the wide end of the homcone is much milder than that of the narrow end, 
oscillating within 23% above or below the value of the reference pressure p 0 . 


The ratio of the amplitudes of the second and the first harmonics of the pressure, p 2 / p s , 
is plotted in Fig. 2 against the reduced amplitude of the first harmonic, pj p 0 , measured at one 

end of the cylinder and at the narrow end of the homcone. The figure shows that the second 
harmonic reaches its maximum at small amplitude of the first harmonic for the consonant 
(cylinder) resonator; however the second harmonic pressure increases more slowly than the first 
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harmonic for the dissonant (homcone) resonator. This indicates that the energy is contained in 
the first harmonic in the homcone resonator, thereby creating a large-amplitude pressure wave in 
the dissonant resonator. In contrast, the energy is efficiently transferred in a cylindrical resonator 
from lower to higher harmonics of pressure at relatively low total pressure amplitudes, thereby 
preventing the formation of large-amplitude pressure waves. 

The amplitude of the first harmonic corresponding to the pressure in the homcone 
resonator is plotted against the frequency of the oscillating resonator in Fig. 3 for two different 

levels of acceleration, ,4 = 2x10^ and A - 1 xlO' 3 . For the smaller acceleration, the pressure is 
uniquely determined at each frequency; for the larger drive amplitude, the pressure takes one of 
the multiple values near the resonance, depending on the direction of the change in frequency. 
The existence of hysteresis and hardening resonance in homcone resonators requires that the 
largest pressure compression ratio be obtained through an upward frequency sweep. The 
numerical simulations show that the frequency increment size near the resonance must be small 
for the convergence of the Multiple Shooting method. 


B. Comparison with experiments 

While the numerical model by itself provides insight into the physics of the nonlinear 
acoustic standing waves in oscillating resonators, experimental verification is required to ensure 
that the model behaves as expected. Experimental data is used to generate the inputs to the 
numerical model. The result of the numerical simulation is then compared to the data gathered 
from the experiment. 
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The acceleration input to the numerical model is generated using measured signals from 
the experimental data. Recall that the model assumes that the resonator is excited using a simple 
sinusoidal acceleration. In addition, the resonator is assumed to be ideally rigid so that the 
acceleration is constant across the length of the resonator body. The arguments made in 

Finkbeiner et al 9 allow the measured acceleration signals to be input to the model by isolating 
the fundamental acceleration harmonic as the driving function for the model. 

The effective viscosity parameter G was originally unknown. The value for G is based in 
part upon the second coefficient of viscosity, due to the high-pressure gradients and frequencies 
present in the acoustics. Previous work 2 assumed a value of G - 0.0 1 , although no justification 
was given for this value. To estimate this parameter, the pressure from the model computed for 
several different values of G, and the results for the varying values of G are then plotted together 
with the experimental data to determine the best match. 10 

In this study, a single comparison between experimental and numerical data is presented 
as an example of the procedure used. The experimental configuration is shown in Fig. 4. A horn- 
cone resonator is rigidly mounted wide-end down on a Labworks ET-127 shaker system. The 
shape of this resonator is given by Eq. (10) with the shape parameters S 0 = 0.028333, 

5, =5.7264, andS 2 =0.25. 

PCB Piezotronics 353B03 accelerometers are mounted at both ends of the resonator. 

PCB Piezotronics dynamic pressure transducers are mounted in the sidewall of the resonator at 

axial positions of X - 0.05 and X = 0.95 . Druck PDCR 130 static pressure transducers and 

thermocouples are mounted at the same axial positions as the dynamic pressure transducers (see 

Fig. 4). The dynamic transducers are optimized for high frequency pressure measurements and 
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are incapable of measuring static pressures due an exponential decay in voltage. The static 
transducers are used to measure this static value and complete the pressure measurement. 

Several acceleration levels are used to excite the resonator; for the purposes of this study, 
an acceleration amplitude of roughly 50 g measured at the narrow end of the resonator is 
presented. At the acoustic resonance of the working fluid, the acceleration signals measured 
near the two ends have small difference in amplitude as graphed in Fig. 5. The dimensionless 
acceleration based on the signal measured at the wide and the narrow ends of the resonator is 

given byA = 9.76xl0' 5 and A = 1.03xl0“ 4 respectively. The difference in the amplitude at the 
two ends is due to elasticity of the resonator. The dynamic pressure frequency response 
measured at the narrow end of the horn-cone resonator is displayed in Fig. 6 along with the 
curves which are computed for varying values of G with a dimensionless acceleration level, 

A = 1 .03 x 1 0" 4 , matching the acceleration measured at the resonator narrow end. The 
computational pressure amplitude is measured at the axial location of X = 0.05 . The result 
shows that the model over-predicts the resonant frequency of the horn-cone resonator by roughly 
2.6%. However, in terms of amplitude and wave form, the model predicts the pressures 
generated in the resonator very well. The matching value of G is between 0.01 1 and 0.012 with 
acceleration referenced at the narrow end. Similar comparison shows, when the computed 
pressure frequency response curves are calculated based on the dimensionless acceleration at the 
narrow end, A = 9.76x1 O' 5 , the matching value of G is found to lie in (0.010, 0.01 1). These 
values of G are, for all intents and purposes, close to those assumed in previous studies 2 . 
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C. Optimal conical shapes 

In following subsections, the results are shown from optimizing each type of resonator 
shapes to achieve maximum compression ratio R c at one end of the resonators. Again, in the 
optimization process the ratio of specific heats is held at y=l .2, the attenuation coefficient 
G=0.01, The acceleration A for a given shape is calculated from Eq. (8), where the value of the 
fixed input power is deduced from that for a reference conical resonator accelerated at the 

amplitude A = 5 x 1(T 4 . The conical resonator contour can be written as 

R(X)=S 0 + S^, for 0 < A' < 1 (11) 

The conical resonator studied in Lawrenson et al and Ilinksii et al. 2 with 50=0.032941 and 
5'i=0.26800 is chosen to be the reference resonator. For the reference resonator, the compression 
ratio reaches its maximum value of R c = 5.0475 when the oscillation frequency is raised to 
fJ = 1.3134 . 

Using the optimization procedure described in Sec. Ill and starting with the reference 
conical shape, the optimal conical resonator is found to have the shape parameters 
S 0 = 0.020833 and 5, = 0. 1 5438 . The compression ratio R c of the conical resonator reaches a 

value of 6.7213, which is about 33% higher than that of the reference resonator. Recall that, due 
to hysteresis, in order to find R c , the entire branch must be traced by incrementing the frequency 
Q. The compression ratio reaches the value when the frequency is increased to 1 .3 123. The 
corresponding acceleration that achieves such a compression ratio has the amplitude 
A = 6.6859 x 1 O’ 4 . According our assumption (8), accelerating the optimal resonator at this 


16 



Li, JASA 


acceleration and frequency consumes the same amount of input power that oscillates the 
reference conical resonator at A = 5 x 1 (T 4 and its resonance frequency Q = 1 .3 1 34 . As stated 
before, to find realistic dimensions of the resonator, the optimization scheme limits the lower 
bound for the dimensionless radius of the resonator at the narrow end to be R( X = 0) = 

0.020833. From the dimensions of the optimal conical resonator, the smaller narrow end, S 0 , 
generates a larger compression ratio. Due to the fixed input-power constraint, the slope of the 
optimal cone S , is finite. The reference resonator and the optimal conical resonator are shown in 

Fig. 7. The second cone (Fig. 7(b)) has the smaller narrow end and smaller slope than the first 
(Fig. 7(a)). 

D. Optimal horn-cone shapes 

The hom-cone geometry used in this study is described by Eq. (10). The three shape 
parameters S 0 , 5, and S 2 are optimized. The first optimization attempt starts with the 

dimensions of the hom-cone given in Lawrenson etal.\ S 0 = 0.028333, 5, = 5.7264 and 
S 2 = 0.25 . Using the original shape parameters, the compression ratio R c = 13.564 is achieved 

when the hom-cone resonator is accelerated at A = 5.06 13 x 1 0~ 4 and the frequency is increased 
to Q = 1 .4674 . In searching for the optimal homcone design, the lower bound of the radius of 
the narrow end S Q is limited to be 0.020833 and the upper bound of S ] to be 6 for practical 

reasons. The separation point between the hom section and the cone section, S 2 , could be any 
value between 0 and 1 . The result of optimization finds that the compression ratio reaches a 
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maximum value of 21.547, 58% higher than that of the homcone of Lawrenson et al\ when 
S 0 = 0.20833, 5, = 5.7263 and S 2 = 0.24893 . This optimal homcone achieves the value of the 

compression ratio when its acceleration has the amplitude A = 5.9782 x 10" 4 and it is oscillating 
at its resonance frequency Q = 1 .4696 . These values of the acceleration and the frequency 
satisfy the requirement of fixed input power formulated in Eq. (8). 

The optimization scheme BFGS is designed for finding a local extreme of a multivariable 
function. Our numerical simulations indicate that the compression ratio, as a function of the 
shape parameters, usually has many local extrema. Our second optimization attempt started with 
a different initial geometry of the homcone and a much higher value of the compression ratio 
R c = 40.203 is obtained when S 0 = 0.020833, S', = 5.5232 and S 2 = 0.22747 . For this 

compression ratio, the resonator is oscillated at the acceleration amplitude A = 6.581 7 xKT 4 and 
the frequency Q = 1 .4259 . Comparing with the compression ratio for the homcone in 
Lawrenson et al } (shown in Fig. 8(a)), the second optimal homcone (shown in Fig. 8(b)) 
improves the compression ratio by more than 196% at the same value of input power. For 
homcone shapes, the principle that higher compression ratio is obtained with smaller narrow end 
is valid, as in the case of conical shapes. Comparing the dimensions of the two homcones, the 
optimal homcone is smaller and the cone section is relatively longer than that of the homcone 
used in Lawrenson et al. 1 . In Fig. 8(c), the pressure waveform at the narrow end is presented for 
the two homcones. The waveform for the new homcone design is more complicated exhibiting 
two notable peaks per cycle. One pressure peak is significantly higher than the other. The two 
pressure waves corresponding to the two hom-cones also have different phases. Under the 
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prescribed conditions, the peak of the pressure is about 4.7 times the ambient pressure in the 
optimal homcone while the peak is about 3.2 times for the homcone in Lawrenson et al. ] . 

E. Optimal cosine shape 

The cosine shape is defined as 

R{X) = S 0 + S, (1 - cos(S 2 *)), for 0 < X < 1 (12) 

Assuming the input power is entirely determined by the interior gas pressure, Chun and Kim 5 
compared three individual resonators of the same volume: one conical shape, one '/ 2 -cosine shape 
(S 2 = n ) and one 3 /4-cosine shape (S 2 =3n/2). They found that the '/ 2 -cosine resonator 
produces the highest compression ratio among the three resonators. Starting with the ‘/ 2 -cosine 
dimensions reported in Chun and Kim 5 , S 0 = 0.025, S l = 0.095, S 2 = n , all three shape parameters 
are optimized using our criteria. The resulting optimal cosine shape parameters are given by 
S 0 = 0.020833, S, = 0.073545 and S 2 = 3.2603 , and the resonator shape is shown in Fig. 9(b). 

For the optimal shape, we obtain the compression ratio R c =13.693 at the frequency 

Q = 1 .4830 , which is about 20% better than that of Chun and Kim 5 at the same level of input 
power. Using different initial guesses of shape parameters for cosine resonators, the 
optimization scheme produced results that yielded lower compression ratios. 

Under the specified conditions (the same power input <W)>, the ratio of specific heats y 
and the viscosity-related parameter G ), the hom-cone shape is found to be better than the cosine 
shape in generating higher compression ratio at the narrow end, opposite the findings of Chun 
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and Kim 5 . As shown in Fig. 9(c), the pressure waveform of the optimal cosine design is similar 
to that of Chun and Kim 5 with a greater pressure peak and a different phase. 

V. CONCLUSIONS 

A local optimization scheme is presented for finding the resonator shapes that maximize 
the pressure compression ratio at one end of an oscillating acoustic resonator. The optimal 
dimensions are reported for cone, homcone, and cosine shaped resonators. The results are 
summarized in Table I, including the shape parameters, corresponding acceleration amplitude, 
resonance frequency, and compression ratio. 

For each type of resonator, a smaller narrow end is found to give a larger pressure peak- 
to-peak ratio. This finding suggests that the number of shape parameters in an optimization can 
be reduced by setting the dimension of the narrow end fixed at a value as small as possible. For 
the types of homcone and cosine shapes, there are many different designs that achieve local 
extrema. Using different initial guesses for the optimal design, the results show that as much as 
1 96% improvement on the compression ratio can be achieved with a fixed level of input power. 
For the shapes considered herein, the homcone shape is found to generate the highest 
compression ratio. The acoustic field in a resonator is a continuous function of the resonator 
shape, the ratio of specific heats and the viscosities. Consequently, the dimensions of optimal 
resonator shape reported herein will undergo small changes if ratio of specific heats and/or the 
viscosities change by a small amount. Searching strategies for globally optimal shapes are under 
investigation. 
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Table I. Summary of the optimization. The shape parameters of the resonators in previous studies and the 
optimized resonators and their compression ratios are listed. The corresponding acceleration levels and 
oscillating frequencies that achieve the compression ratios are also shown. 



So 


S 2 

A (xlO' 4 ) 

Q 

K 

Cone: 







Ilinskii, ei al. 

0.032941 

0.268 

N/A 

5.0 

1.3134 

5.0475 

Optimal 

0.020833 

0.15438 

N/A 

6.6859 

1.3123 

6.7213 

Horn-cone: 







Lawrenson, et al. 

0.028333 

5.7264 

0.25 

5.0613 

1.4674 

13.564 

Optimal 

0.020833 

5.5232 

0.22747 

6.5817 

1.4259 

40.203 

Cosine-shape: 







Chun and Kim 

0.025 

0.095 

3.1416 

5.6299 

1.4837 

11.440 

Optimal 

0.020833 

0.073545 

3.2603 

5.0613 

1.4830 

13.693 
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Figure Captions 

Figure 1 . The reduced pressure p / p 0 at the ends of a homcone and a cylinder at their 
corresponding resonance frequencies for the same acceleration A = lx 1(T 3 . 

Figure 2. The ratio between the amplitudes of the second and the first harmonics of the pressure 
wave, p 2 ! P\, is plotted as a function of the dimensionless amplitude of the first harmonic 
reduced pressure pj p 0 . These amplitudes are calculated for the narrow end of the homcone 
and at one end of the cylinder. 

Figure 3. The effect of the frequency and the amplitude of acceleration on the dimensionless 
amplitude of the first harmonic reduced pressure /?, / p 0 . The pressure is calculated for the 

narrow end of the homcone for the acceleration A = 2x10^ (crosses) and A = 1 x 10' 3 (circles and 
squares). 

Figure 4. Experimental apparatus. 

Figure 5. The dimensionless acceleration signals measured at the two ends of the homcone 
resonator by PCB accelerometers. 

Figure 6. The dynamic pressure frequency response measured at the narrow end of the hom-cone 
resonator and the pressure responses predicted by the numerical model computed for varying 
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values of G with a dimensionless acceleration level, A = 1 .03 x 1 0 4 , matching the acceleration 
measured at the resonator narrow end. 

Figure 7. Conical shapes: (a) the reference conical resonator; (b) the optimized conical 
resonator. For the same input power, the compression ratios corresponding to the shapes in (a) 
and (b) reach the values of 5.0475 and 6.7204 respectively. For shape parameters, see the text. 

Figure 8. (a) The homcone shape in Lawrenson et al. ] (b) Optimized homcone shape. For shape 
parameters, see the text, (c) The pressure waveform at the narrow end for the homcone in 
Lawrenson et al . 1 (the solid line) and that for the second optimized homcone (the dashed line) 
are shown. 

Figure 9. (a) The A-cosine shape as in Chun and Kim 5 , (b) The optimized A-cosine shape, (c) 
The pressure waveform at the narrow end for the A-cosine resonator in Chun and Kim 5 (the solid 
line) and that for the optimized A-cosine shape (the dashed line) are shown. 
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